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ABSTRACT 

The Extended Narrow Line Region (ENLR) of Active Galactic Nuclei (AGN) is 
a region of highly ionized gas with a size of few up to 15-20 kpc. When it shows 
a conical or bi-conical shape with the apexes pointing towards the active nucleus, 
this region is also called ionization cones. The ionization cones are an evidence of 
the Unified Model that predicts an anisotropic escape of ionizing photons from the 
nucleus confined to a cone by a dusty torus. Many details about the complex struc- 
ture of the ENLR still remain unveiled, as for example the origin of the ionized gas. 
Here we present new results of a study of the physical and kinematical properties of 
the circumnuclear gas in the nearby Seyfert 2 galaxy NGC 7212. Medium and high 
resolution integral field spectra and broad-band photometric data were collected and 
analysed in the frame of an observational campaign of nearby Seyfert galaxies, aiming 
to handle the complicated issue of the origin of the gas in the ENLR. This work is 
based on: (i) analysis of gas physical properties (density, temperature and metallicity), 
(ii) analysis of emission line ratios, and (iii) study of kinematics of gas and stars. By 
reconstructing the [Olll]/H/3 ionization map, we pointed out for the first time the 
presence of an ionization cone extended up to about 6 kpc, made by a large amount of 
low metallicity gas, kinematically disturbed and decoupled from stars, whose highly 
ionized component shows radial motions at multiple velocities proved by the complex 
profiles of the specral lines. Since NGC 7212 is a strongly interacting triple galaxy 
system, the gravitational effects are likely to be at the origin of the ENLR in this 
Seyfert galaxy. 

Key words: galaxies: Seyfert - galaxies: individual: NGC 7212 - techniques: spec- 
troscopic - line: profiles 



1 INTRODUCTION 

Ionizati on cones are an ev idence of the validity of the Unified 
Model (lAntonuccilll993 h. which postulates an anisotropic 
escape of photons from the Active Galactic Nuclei (AGN) 
confined to a cone by a dusty torus. Up to now, a bright 
Extended Narrow Line Region (ENLR) was found in only 
25 galaxies. Many details about the complex structure of the 
ENLR still remain unveiled. The small number of ionization 

~^ This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile, and it is 
partly based on observations made with the 6-m telescope of the 
SAO RAS operated under the financial support of the Ministry of 
Science of the Russian Federation (Registration Number 01-43). 
t E-mail: valentina.cracco@unipd.it 



cones could be related to the still open problem of the origin 
of the ionized gas in the ENLR. In fact, the ENLR gas could 
be part of the interstellar medium (ISM) of the host galaxy, 
photoionized by the active nucleus, or material ejected by 
the nucleus in strong interaction with a radio-jet. Otherwise, 
the gas could be acquired from the intergalactic medium, if 
the galaxy is in a dense environment or by means of grav- 
itational interactions, in case of minor merger events. This 
issue is strongly related to the mechanisms of AGN feeding. 

Within this frame, we are carrying out an observational 
program taking advantage of the integral-field spectroscopic 
technique to investigate the physical and kinematical prop- 
erties of the ENLR in few nearby Seyfert 2 galaxies (z < 
0.03) with possible ionization cones or at least extended 
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[O III] emission. Here we present new results of the first stud- 
ied object, NGC 7212. 

NGC 7212 (z=0.0266) is a Seyfe rt 2 galaxy in a n 
interacting system of three galaxies (|Wasilewskil Il98l[ ). 
iDurret WarinI (|l990l l by means of optical longslit spec- 
troscopy in the H/5 + [O ill] spectral range, found ionized 
gas extending for ^ 17 arcsec along PA = 37° and 127°, and 
a high excitation value R=I([Oiii])A 5007+4959/I(H/3) = 19 
in th e nuc l eus, ranging from 5 to 28 in the ionized emis- 
sion. iTranI (| 19951 ) detected a jet-like highly ionized feature 
extended up to 10 arcsec from the nucleus at PA = 170° in 
ground-based [O III] and Ha images. This feature is exactly 
parallel to the axis of the small s cale double radio source (0.7 
arcsec separation) published bv lFalcke, Wilson, &; SimpsonI 
(| 19981 ), who did not find a clear evidence of an ionization 
cone, since the Northern side of the galaxy is partially ob- 
scured by dust lanes. However, they suggested that the dis- 
tribution of excitation to wards the South is consistent with 
a ragged ionization cone. iKotilainenI (| 19981 ) published the 
B — I map of NGC 7212 showing a very blue and fan-shaped 
emission region extending from the nucleus towards South 
(PA = 165°) with total size of 2.3 arcsec, and a dust lane 
situated on the other side of the nucleus, at 3.7 arcsec and 
PA 280°. Hubble S pace Telescope (HST) [O iii] images by 
ISchmitt et aD (|2003l ) showed that the [O III] emission is ex- 
tended up to ~3 arcsec from the nucleus along PA = 170° 
with dimensions 2.1 x 4.8 arcsec. The emission is diffuse and 
composed of sever al individual kno t s Nor thwards and South- 
wards the nucleus. iRaimann et al.l (|2003h found that the nu- 
clear stellar component of NGC 7212 is dominated by 10- 
Gyr metal rich (solar or above solar metallicity) stellar pop- 
ulation, with a contribution of 15 per cent of the total flux 
from the 1-Gyr component and 15 per cent f rom the 3-Myr 
compo nent or from a featureless continuum. iBennert et al.l 
(|2006l ) confirmed the [O ill] emission along PA = 170°, but 
extended up to 12 arcsec from the nucleus, i.e. four times 
larger than the extension seen in the HST image in the same 
direction, but smaller t han the maximum extent reported by 
iDurret Warin (1990"), who observed at different position 
angles. The reddening corrected excitation value in the cen- 
tral spectrum is R ~16 and varies between 6 and 17 in the 
central 24-arcsec-wide region. The reddening in the centre 
is rather low (E(B— V)= 0.33 ± 0.01 mag) and decreases to 
a value of ^0.07 mag at 1 arcsec North- Westerly of the nu- 
cleus. On both sides of this region, it increases and reaches 
its maximum value at 4 arcsec South-Easterly and 7 arc- 
sec North- Westerly of the photometric centre (AE(B— V)^l 
mag) . These maxima ma y be attributed to dust lanes seen i n 
the continuum image bv lFalcke, Wilson, &: SimpsonI (|l998l ) . 



2 OBSERVATIONS AND DATA REDUCTION 

NGC 7212 was observed at the Russian 6-m telescope of 
the Special Astrophysical Observatory (SAO RAS) with 
the Multi Pupil Fiber Spectr ograph (hereafter MPFS, 
lAfanasiev, Dodonov, MoiseevI [2001). This integral field 
unit takes simultaneous spectra from 256 spatial elements, 
the field of view (FoV) is 16 square arcsec with a spatial 
sampling of 1 arcsec. Spectra of NGC 7212 were obtained in 
September 2005 in the range 3800-7300 A with a 600 lines 
mm~^ grating, and in July 2007 in the range 4000-7300 A 




Figure 1. i?-band image of NGC 7212 obtained at the SAO 
6-m telescope with SCORPIO. The white (black) square is the 
MPFS FoV for the low (high) resolution data. The FoV is 16x16 
square arcsec, corresponding to 8.2x8.2 kpc^ (scale = 0.513 kpc 
arcsec" Hq = 75 km s"-*^ Mpc"-*^). North is up and East is on 
the left. 



Table 1. The observational and archival data 



Date Instrumental Spectral range AA 5\ Tex] 

(A) (A) (A px-i) (s) 



Sep 2005 


MPFS 


3780-5800 


6 


1.46 


4800 


Sep 2005 


MPFS 


4340-7310 


6 


1.47 


5400 


Jul 2007 


MPFS 


4140-5650 


3 


0.76 


4800 


Jul 2007 


MPFS 


5765-7260 


3 


0.76 


1800 


May 2010 


MagE 


3600-7000 


0.85 




7200 


Sep 2008 


SCORPIO 


B 






600 


Sep 2008 


SCORPIO 


V 






660 


Sep 2008 


SCORPIO 


R 






510 


Sep 1994 


WFPC2 


F606W 






500 



with a 1200 lines mm~^ grating. The data have dispersion of 
1.5 A px~^ and 0.8 A px~^, and instrumental resolution of 
^-^6 A and ^3 A, respectively. NGC 7212 was also observed 
with the high resolution (R^8000) MagE echelle spectro- 
graph at the Magellan Telescopes of the Las Campanas Ob- 
servatory (Chile), with the slit oriented along the ionization 
cones (PA-170°). From HST-MAST archive a WFPC2 im- 
age in the F606W broad band filter was extracted. Broad 
band BVR images were obtained in Septe mber 2008 at the 
SAO 6-m telescope with SCORPIO ( Afanasiev Mo iseevI 
l2005h (Fig. [1]). All the observational and archival data are 
listed in Table [1] where AA is the spectral resolution in A 
and 6X is the dispersion in A px~^. 

The MPFS data were reduced using P3d ( Sandin et al.l 

l2010h which is a general data-reduction software written in 
IDL language and developed to work with fibre-fed integral- 
field units (IFU) of any integral-field spectrograph. Up to 
now, it has been configured and tested with PMAS, VIRUS- 
P, SPIRAL and MPFS. P3d makes the reduction steps up 
to wavelength calibration automatically, and allows to in- 
teractively inspect and optimize parameters when required. 
It also provides graphical tools to check the raw data and 
the output of the different tasks. In order to perform the 
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flux calibration, we observed a spectrophotometric standard 
star for each run. Its spectra were extracted, flat fielded and 
wavelength calibrated with P3d, then sky subtracted and 
summed together, to collect the total flux, with iraf. The 
usual IRAF tasks for the flux calibration were applied to the 
data. The echelle spectrum was reduced with iraf following 
the standard procedures (bias subtraction, flat-field correc- 
tion and wavelength calibration). No other calibrations or 
corrections were applied. Before reducing this spectrum, the 
orders were extracted using APALL with the option strip. 

After the reduction procedures, the IFU spectra were 
corrected for Galactic reddening, using the iraf task dered- 
DEN and the value of F-band absorption (Ay = 0.238) 
derived from NED (NASA/IPAC Extragalactic Database), 
then we corrected for telluric absorptions. In fact, at z = 
0.0266 the [Sii] doublet falls in the range 6880-6930 A, and 
therefore it is only partially absorbed by the atmospheric B 
band (6860-6890 A). We used two iraf scripts, ATMO and 
RMAT, written by our group. ATMO fits the continuum of the 
flux- calibrated standard star and divides the fitting by the 
stellar spectrum. Then, all the intensity values at a A shorter 
than 6700 A are put equal to 1. Finally, the corrected spec- 
trum is obtained using the task rmat which simply mul- 
tiplies the observed spectrum by the output of ATMO. We 
also applied a correction for atmospheric refraction. We used 
two iraf scripts written by us, named ACORR and darc- 
SPEC. The first task determines the correction to apply for 
atmospheric refraction, the second one applies this correc- 
tion. In particular, ACORR creates images of the emission for 
20-30 intervals centred at different A along the dispersion 
direction, then it calculates the centroid of each image using 
IMCENTROID and allows to interactively fit with a straight 
line the relations A vs. x and A vs. using nfitId. darc- 
SPEC uses the x{X) and y{X) functions and for each A along 
all the spectral range, makes an image, applies the needed 
shift (taking as a reference the centroid at A = 5500 A) and 
then transforms the image in a table with three columns: x, 
2/, A. From these output tables the spectrum corrected for 
atmospheric refraction is reconstructed. 

In order to accurately measure the emission lines, in 
particular hydrogen and helium Balmer lines, the underly- 
ing stellar component was subtr acted from the spectra. We 
appl ied the starlight software (|Cid Fernandes et al ] l2007l . 
I2OO5) to fit the galactic stellar component. Before being 
analysed with starlight, the spectra were shifted to rest- 
frame wavelengths with newredshift and their dispersion 
was modified to a value of 1 A px~^ with dispcor. We used 
45 synthetic spectra by combining 15 ages (from 10^ yr up 
to 13 X 10^ yr ) with 3 metalhcities (Z 0.004, 0. 02 and 0.05 
Zo), and the lCardelh. Clavton. Mathid (| 19891 ) (hereafter 
COM) extinction function. The synthetic spectrum of the 
stellar contribution obtained from starlight for each fibre 
was subtracted from the observed one to obtain a pure emis- 
sion line spectrum. These procedures were performed only 
for low resolution spectra, because high resolution data were 
used only to study the kinematics. Moreover, in high resolu- 
tion spectra the Balmer absorption lines could be neglected 
because of the low S/N ratio (less than 5) of the continuum. 

For the analysis of IFU data, we used pan (Peak ANaly- 
sis) , an IDL software that fits emission lines with a graphical 
interface, based on Craig Markwardts mpfit. We applied a 
single Gaussian function for all the emission lines except Ha 




Figure 2. Left: the S/N of the IFU low resolution data mea- 
sured at 5500 A (solid line) and at 7000 A (dashed line). Right: 
the histogram of the x^, evaluating the quality of the starlight 
fitting. 

+ [N 11] and [S 11] for which we used custom functions which 
allowed to fit all the lines simultaneously, pan is useful for 
integral field data since it can read multiple spectra, and fit 
the initial guess to all the spectra automatically. The spec- 
trum, fitting and residuals can be reviewed, being displayed 
on-screen. The initial parameters (amplitude, position and 
width) can be specified interactively and constrained to en- 
sure a reasonable result: it is possible to fix to the same value 
the FWHM of the doublets, like [N 11] and [S 11] . For the high 
resolution data of NGC 7212 we used the midas package 
XAlice because it was necessary to fit the [O ill] emission 
lines with two components, a broad and a narrow one. 

The reliability of our measurements was verified by cal- 
culating the signal to noise ratio (S/N) of the continuum 
for MPFS data (see Fig.[2])- These data have small values of 
S/N: the median values for low resolution data are ~5.6 and 
^^5.4, measured at rest-frame 5500 A and 7000 A. The 
evaluates the quality of the starlight fitting: it is peaked 
around 1.6 (see Fig. [2}, instead of 1 probably because the 
S/N ratio of the analysed spectra is not sufficiently high. 
We calculated the flux errors, assuming that the determina- 
tion of the position of the continuum level gives the main 
contribution to the errors when measuring the emission line 
fluxes. 



3 PHYSICAL PROPERTIES OF GAS AND 
STARS 

All the measured fluxes were corrected for the internal red- 
dening, using the Balmer decrement and assuming the theo- 
retical value of the intensity ratio between Ha and H/3 equal 
to 2.86. The visual absorption, A(y), was estimated by ap- 
plying the CCM extinction law. Once we obtained A(V), 
we corrected all the observed fluxes. In order to study the 
spatial distribution of the extinction, we reconstructed the 
A(V) map (see Fig. [3]), which is clearly limited to the re- 
gion where the H/3 emission line can be measured. In the 
nucleus the reddening is low, A(l/) = 0.78 mag, higher val- 
ues, A(y) = 1.6-2 mag, are found at about 7 arcsec North- 
West and at about 4 arcsec East of the nucleus, in agree - 
ment with the reddening obtained bv lBennert et al.1 (I2OO6I), 
E(B-V) - 0.33 mag in the centre and E(5 - F) ^ 0.8 mag 
in the northwestern and southeastern regions (see their fig. 
3, upper-right panel). We estimated the extinction from the 
stellar spectra by means of starlight, which also gives as 
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Figure 3. Left: A(V) map, obtained by applying the CCM ex- 
tinction law to the ratio of the intensity of Ha and H/3, with the 
stellar continuum contours overlaid (in black). Right: The his- 
togram of the A(V) values obtained for gas (solid line) and stars 
(dashed line). 




Figure 4. Left: the [Olll] emission line map with the overlaid 
contours of the stellar continuum emission at 5500 A in black. 
Right: the same for Ho;. The gas emission is oriented along PA = 
0° and the stellar continuum is oriented along PA = 45°. North 
is up and East is on the left. 



output the A{V) values for each spectrum. By comparing 
the A{V) obtained from gas and stars, we can see that the 
distributions are peaked at different values, larger for the 
gas: the median value for gas is 1.38 mag and for stars is 
0.68 mag. If the ENLR gas were dust-free, then the A{V) 
measured by stars and gas would be the same. The larger ex- 
tinction observed in the gas component suggests that there 
is also a fraction of dust associated to the ENLR. 

Up to 30 emission lines were detected and measured at 
least in the central regions, within a radius of about 4 arcsec. 
The maps of the brightest emission lines ([O ii] A3727, [O ill] 
A4959, [Oiii] A5007, [Oi] A6300, [Nil] A6548, Ha A6563, 
[Nil] A6584, [Sii] A6716, [Sii] A6731) show an elongated 
shape: the ionization gas is extended up to 15 arcsec (cor- 
responding to ^7.7 kpc, scale = 0.513 kp c arcsec"^) with 
a PA :^ 0°. A similar size was detected bv lDurret &: WarinI 
(Il990h with longslit spectroscopy. They found ionized gas 
extending for about 17 arcsec at PA = 37° and PA = 127°, 
corresponding to the minor and major axis of the galaxy, re- 
spectively. We verifie d that this emission is within our FoV. 
iBennert et alJ (|2006l ) found ionized gas extending out to 12 
arcsec along a PA ^ 170°, as well. We compared our emission 
lin es maps of [ O ill], Ha and [Nil] with the images published 
bv lTranI (|l995l ). Fig.[5]shows a good match between the con- 
tours of these images and our two-dimensional (2D) maps. 
The emission is not oriented as the stellar continuum emis- 
sion which has the major axis at PA = 45°, in Fig. |4]the 
[O III] and Ha emission line maps are shown with the con- 
tours of the stellar continuum emission at 5500 A overlaid. 
[O III] A4363 and He ii A4686 maps are clearly elongated as 
well with PA ^ 0° and extended up to 9-10 arcsec (corre- 
sponding to ^ 4.5-5 kpc). Fainter emission lines are visible 
not only in the very inner parts, where the influence of the 
ionization source is more important, but also at 3-6 arcsec 
(1.5-3 kpc) from the nucleus. The coronal emission lines of 
iron with different ionization degree were detectable as well 
and were measured, but only in the central spectra, probably 
because here both the ionization and the S/N are sufficiently 
high. We detected [Feiii] A4658, [Fev] A4228, [Fevii] A5158 
and [Fevii] A5721 and [Fevii] A6087 (Fig. [6]). 

We calculated the density with the temden iraf 
task, by using both the [S ii] A6716/A6731 and the [Ariv] 
A4711/A4740 ratios. The [Sii] and [Ariv] lines have differ- 
ent values of critical density (logA^e = 3.2, 3.6, 4.4 and 5.6, 





Figures. Left: 2D map of the flux distributions for [Olll] with 
the iTranI (Il995l ) contours overlaid in black. Right: the same for 
Ha. North is up and East is on the left. 
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Figure 6. The nuclear spectrum of NGC 7212. In the upper panel 
the brightest emission lines are clearly visible, in the bottom panel 
the Argon and Iron emission lines are shown. 
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Figure 7. Left: histogram of density obtained using a value of 
temperature Te = 10^ K, the median value is 450 cm~^ (dashed 
line). Right: histograms of temperature obtained with [Olll] (solid 
line) and [Sll] (dashed line) ratio, the median values are 1.8 x 10^ 
K and 1.5 x 10^ K, respectively. 



respectively), thus we can assume that these lines are emit- 
ted by gas in different physical conditions. The [Ar iv] ratio 
gives information about higher density and higher ionization 
gas, while lower density gas can be studied by means of the 
[S ii] ratio. We used an input value of temperature T — 10^ 
K and we calculated the density for each spectrum in which 
we measured these ratios. Unfortunately, the [Ariv] lines 
could be measured for few central spectra only, giving not 
significant information about the spatial distribution of the 
high density gas. Therefore we calculated the density using 
the [S li] doublet, which is detectable in a more extended 
region, finding a median value of ^ 450 cm~^ (see Fig. (T)). 
The central density obtained with [Ar iv] is ~ 2.2 x 10^ cm~^ 
(assuming Te = 10^ K) and is in good agreement with the 
value obtained with [S ii] (^ 1.4 x 10^ cm~^). These higher 
values of density ( 1-1. 4 x 10^ cm~^) are in a n internal region 
oriented roughly East- West. iBennert et aP (2006) also esti- 
mated high density values (> lO^cm"^) within an aperture 
of 2 arcsec on the nucleus (see their fig. 6, upper-right panel). 
The gas temperature was obtained by calculating the ratio 
between the flux of the [O ill] AA4959+5007 and the flux of 
the [Oiii] A4363, and by calculating the ratio between the 
flux of [S II] AA6717+6731 and the flux of [S ii] AA4068+4076. 
We determined the temperature with these two ratios, in 
case that the A4363 and the AA4068+4076 emission lines 
could be measured. We derived the temperature using tem- 
DEN and the values of density calculated with the [S ii] ratio. 
The histograms for the two determinations of temperature 
are shown in Fig. [71 The median values are 1.8x10^ K and 
1.5 xlO'^ K for [Oiii] and [Sii] ratios respectively. For high 
ionization gas, we should calculate the density using high 
ionization doublets, but in this case, we detected the [Ar iv] 
emission lines in the inner parts only, and we could estimate 
the density using this ratio only for few nuclear regions, thus 
we used the [S ii] density determination for the temperature 
estimate with the [O ill] ratio. Higher values (up to 2-3 x lO'^ 
K) were found in the northern part of the FoV, while in the 
central regions we measured lower values (1-1.5 xlO^ K) 
with [Oiii], and even less with [Sii] (4-9 xlO^K). 

In case of ionization by a power-law spectrum, as for an 
AGN, the determination of the metallicity is complicated be- 
cause the ionization structure can be complex. There are not 
direct methods to determine the metallicity, and compari- 
son with photoionization models must be used. To estimate 



The origin of gas in the ENLR 5 



■ 1.0 

" □ 1.5 
2.0 

' o 2.5 
3.0 

A 3.5 


Zo 
Zo 
Zo 
Zo 


□ □ mO(5 


1 

o O 

o 










- 






o 






' ■ ■ * ■ ""^ ' ° 
















1 


1 



-0.5 0.5 

Log ([N II]/[S 11]) 



Figure 8. The graph of the log([N ll]/[S ll]) versus 
log([N ll]/[0 ll]) compared with theoretical models by IVaonal 
((20101) with D/G=l. The observed values are black crosses. 

the gas metallicity, we compare t he measured fluxes with 
CLOUDY models from IVaoiiar(|2010l ). using the emission line 
ratios [Nil] A6584/[Oii] A3727 vs. [Nil] A6584/[Sii] A6724. 
Vaona showed that this is the more sensitive to metallicity 
diagram. The models have a low resolution in metallicity 
because he used few Z/Z© values, nevertheless the trend is 
clear. We used models with dust-to-gas ratio (D/G) equal 
to 1. The plot and the comparison between the observed 
values and the theoretical models allowed us to estimate a 
sub-solar metallicity for the ionized gas in NGC 7212 (see 
Fig. ID). 

From the STARLIGHT output, we reconstructed the 2D 
maps of the stellar velocity field, as well as the 2D maps of 
the distribution of the stellar population with different ages 
and metallicities. The maps were built using only the models 
in which the adev parameter was less than 50 per cent. This 
parameter is an output of starlight and corresponds to the 
average percentage value (|Oa — Mx\/Ox) of the deviation 
between the observed spectrum (Ox) and the modelled one 
(Ma) over all fitted pixels. Then we selected only the spectra 
with S/N > 3. We obtained the flux distributions for stellar 
populations in different ranges of age: young (t < 10^ yr), 
medium age (lO^yr < t < 10^-^yr), and old {t > lO^'^yr); 
and with different metallicities: sub-solar (Z = 0.004), solar 
(Z = 0.02) and super-solar (Z = 0.05). We created the B - R 
colour image from the broad-band images obtained at SAO 
and we compared it with the 2D maps showing the stellar 
population distribution. The colour image (Fig. [9l bottom- 
right panel) shows a red ring-like structure North and North- 
West of the nucleus, corresponding to the large dust lane 
already visible in the broad-band images and in our A{V) 
map, and a blue region Sout h-East of th e nucle us. Both these 
structures were detected bv iKotilainenI (| 19981 ), who found a 
blue fan-shaped region along PA = 165° extended about 2.3 
arcsec. Our region has the same size but its shape is circular 
due to our seeing of about 1.5 arcsec, worst than the 0.7- 
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Figure 10. Maps of diagnostic diagrams with the i^-band SAO contours overlaid (in black). Top: [Olll]/H/3 (left) and [Oi]/Hq; (right). 
Bottom: [Nil] /Ha (left) and [Sll]/Ha (right). 



arcsec seeing in the Kotilainen images. As in lRaimann et aP 
(|2003l ) we found that the old population (t > 10^-^yr) is 
the main component and in the nucleus we found also a 
younger component (lO^yr < t < 10^"^yr). In particular, in 
a large ring-like structure the 70-80 per cent of the total light 
is due to the old stellar population. A lower precentage of 
the light (50-60 per cent) is observed in a central elongated 
region oriented at PA = 90°, where the intermediate age 
population accounts for 30-40 per cent of the total. This 
apparent hollow in the distribution of the light of the old 
stars is simply an effect of the presence of younger stars 
with smaller M/L ratio. Finally, about 10 per cent of the 



light seems to be emitted by young population, but we stress 
that it could be an effect of diffuse light from the cone and 
scattered light from the AGN. In fact starlight does not 
distinguish between OB stars contribution and a power-law 
continuum. Therefore, even if in a Seyfert 2 galaxy the AGN 
continuum is strongly absorbed with respect to Seyfert 1, the 
detection of the young and very young population could be 
wrong, due to the contamination from the AGN featureless 
continuum. In addition, high values of intermediate stellar 
population are found in the interaction regions (North-East 
of the nucleus). This component could result from the on- 
going merger. 
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Figure 9. Top: 2D maps of the old (left) and medium age (right) 
stellar population. Bottom: the 2D map of young population (left) 
and the B — R color (right). The cross is the position of the galaxy 
nucleus. North is up and East is on the left. 
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Figure 11. The VO diagnostic diagrams. The squares are regions 
inside the cone, the triangles are regions outside. 



We used the diag: nQstic ratios, applied in the BPT 
and VO di agrams ( Bal dwin. Phillips, TerlevichI Il98ll : 
IVeilleux fc"Osterbrock ..1987f )7 to investigate the ionization 
mechanisms in act in the ENLR. The 2D maps of the diag- 
nostic ratios are shown in Fig. [lOl The [Nil] /Ha, [Oi]/Ha 
and [Sii]/Ha maps were created using the non-reddening 
corrected values, in order to study a more extended field, not 
being limited by the small size of the H/3 emission. From the 
2D maps, we detected for the first time an extended struc- 
ture with high values of [O iii]/H/3, larger than 10, suggesting 
the possible presence of an ionization cone oriented as the 
galaxy minor axis (see Section [5])- The [Oiii]/H/3 ionization 
map shows high values even far from the nucleus. If pho- 
toionization is the only ionization mechanism, the photon 
flux is diluted with increasing distance, and it is not expected 
to remain high at large distance from the nuclear source. We 
speculate that this could be a density effect, in fact the ion- 



ization parameter U 



where Qion is the number of 



ionizing photons per second emitted by the source, depends 
both on the inverse of the density and of the distance. So 
high values of U could be due to high values of ionizing pho- 
ton flux or to low density. Inside the cone, we found low 
density, but without a gradient, the value is approximately 
the same along the whole size of the cone. Thus, we can ex- 
clude that high values of U are due to low density. We can 
invoke the presence of shocks as an additional mechanism of 
ionization. Indeed, some regions inside the cone show elec- 
tron temperature values larger than 3 x 10^ K, suggesting 
that shocks could be at work. However, this hypothesis will 
be tested with photoionization+shocks models (Contini et 
al. 2011, in preparation). 

In addition, according to the very high values of 
[Oiii]/H/3, we can also exclude a contribution from star 



formation. The diagnostic ratios of [Nil] /Ha and [Sii]/Ha 
are not smoothly distributed, showing higher values (respec- 
tively ^1-1.3 and ^0.7-0.9) orthogonal to the cone and out- 
side it. Furthermore, these ratios show higher values in the 
North-Eastern regions, namely in the interacting regions be- 
tween the two galaxies, as clearly seen in Fig. [10] with the 
galaxy contours overlaid. Here the ratios can be strongly in- 
fluenced by interaction effects, such as shocks. The [O i]/Ha 
map is peculiar. It shows an elongated shape along the 
galaxy major axis, with high values oriented as the high 
density structure. High values of [O i]/Ha are related to col- 
lisions. The [O i] emission line is formed in the recombination 
region, where O^, H^ and H exist. This region is more ex- 
tended if the ionizing source has a power-law spectrum. This 
ratio is very useful to distinguish between photoionization 
by AGN or stars. We plotted the VO diagnostic diagrams, 
using the reddening corrected fluxes, both for the regions in- 
side and outside the ionization cones (see Fig. [TT]) . In order 
to separate the different regions in the diagnostic diagram 
(AG N, Hii refflons and LINERs), we used the expressions 
from iKewlev et aP (2006). By analysing these diagrams we 
noticed that all the regions occupy the AGN area and are 
distributed vertically towards the shock region. This means 
that for each region, inside and outside the cone, the ionizing 
source has a power-law spectrum. The only difference is that 
for the regions outside the cone the ionization decreases with 
increasing radius, instead, inside the cone, it remains high 
also far from the nucleus (see Fig. [20]). Possible explanations 
are that the ionizing photons from the cone are diffused, or 
that the cone aperture angle is larger. 
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Figure 12. Left: 2D maps of flux (top), velocity (middle) and 
FWHM (bottom) for high ionization gas ([Olll]). Right: the same 
for low ionization gas (Ha). 

4 KINEMATICS 

For the main emission lines, we also built the 2D maps of 
velocity and FWHM (see in Fig. [12] two examples for [Oiii] 
and Ha) . The measured velocities were corrected for the sys- 
temic velocity, considering the value obtained from the nu- 
clear spectrum, that is the spectrum with the highest flux 
in the continuum map. For each spectrum the mean instru- 
mental FWHM, obtained fitting few lines of the compari- 
son lamp, was subtracted quadrat ically from the measured 
values, in order to take into account differences of the in- 
strumental resolution in the field. The maps of the velocity 
field show values ranging from —180 km s~^ to +150 km 
s"^ for [Oiii] and from -220 km s"^ to +180 km s"^ for 
Ka. The maps of the FWHM, after the correction for in- 
strumental width, show values between 200 km s~^ and 800 
km s~^ for both [O ill] and Ha lines. The velocity maps 
show distortions East of the nucleus, especially in case of 
Ha, [Nil] A6548,6584 and [Sii] A6717,6731 lines, while the 
[O III] A5007 line shows a more regular pattern. In particu- 
lar, the Ha velocity map in the eastern region has velocity 
values between 130 km s~^ and 200 km s~^. This zone ap- 
pears to be located towards the interaction region. Its high 
FWHM values (300-600 km s"^) could be an effect of the 




Figure 13. The velocity maps of [O ill] (left), Ho; (middle) and 
stars (right) with the continuum contours overlaid. 



gravitational interaction between NGC 7212 and its com- 
panion galaxy. However, in the HST image and in the SAO 
broad-band image and color, this region shows no particular 
structure and has a low surface brightness (10~^^-10~^^ erg 
cm-2 s-^). In the [Oiii], [Oi], [Nil], Ha and [S ii] FWHM 
maps it is clearly visible an elongated structure, with higher 
values, up to 700-800 km s~^, oriented at about 90° with 
respect to the maj or axis of the emission . The FWHM of 
[Oiii] published by Durret k Warin (1990) ranges from 320 
to 385 km s~^. This disagreement is likely due to the orienta- 
tions of their slit which do not correspond to this elongated 
structure. In any case, we observed values between 500 and 
600 km s~^ in the nuclear region. This feature has been al- 
ready observed in other Seyfert galaxies, like Mrk 3 (Di Millj 
l2007l ) Mrk 34, Mrk 1066, Mr k 348, Mrk 1, NGC 2992, NGC 
5728 (jStoklasova et al.ll2009l ). The explanation of this effect 
is still under debate. This structure is also found when we 
fit the high resolution [O ill] emission line profile with two 
separate components. The minor axis of the stellar veloc- 
ity field is aligned with the photometric minor axis, on the 
contrary the gas velocity field is rotated with respect to the 
stellar one (Fig. I13p . From the analysis of the low resolu- 
tion data we can see that stars and gas are not aligned, the 
[O III] emission is elongated along a direction orthogonal to 
the major axis of the galaxy. This could imply that the gas 
is not co-planar with the stars. It is interestin g; to comment 
that the jet-like structure, studied bv iTranI (|1995l ). extend- 
ing up to about 10 arcsec and oriented at PA ~ —10° and 
showing redshifted velocity of 180 km s~^, is not consistent 
with material ejected by the nuclear source, but instead it is 
rotating ionized gas as clearly visible from our [O III] veloc- 
ity map. The density map is aligned as the FWHM maps. 
By comparing the high density (A^e = 700-1400 cm~^) re- 
gions in the map to the FWHM maps (see Fig. [T4|) . we found 
high values of FWHM (- 500 km s"^) in regions showing 
high density. These are hints of compressed gas, or stream- 
ing of gas, and can be consistent with shocks. The North- 
Eastern region with higher FWHM and a peculiarity in the 
velocity field is nearby the region of the merge between the 
galaxies and not connected to the ionization cone visible in 
the [Oiii]/H/3 map. All the previous peculiarities could be 
caused by this interaction (see Fig. [T4|) . From the analysis 
of the emission line profiles it is possible to point out the 
presence of asymmetries that can be due to additional and 
non-rotational kinematical components. With the low reso- 
lution data, we were not able to identify and analyse mul- 
tiple components in the spectra line profiles, although some 
asymmetries in the line profile for the higher S/N spectra 
were observed and, as a consequence, the fitting was more 
difficult to be performed. From the higher resolution IFU 
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Figure 14. The maps of electron density with the Yia FWHM 
contours overlaid (left) and of Ho; FWHM with continuum con- 
tours overlaid (right). 



data, initially we fitted [O ill] and Ha emission lines with a 
single component and we built the 2D maps of flux, veloc- 
ity and FW HM. We used these maps as input for rotcur 
(|Begeman| [l989). the task of the GIPSY (Groningen Image 
Processing System) software, to model the 2D map and to 
obtain the deprojected values of the velocity in order to build 
the rotation curve. This task uses tilted-ring models: it mod- 
els the velocity map using concentric rings and taking into 
account the inclination on the line of sight. Before starting, 
we needed to calculate the inclination angle (i) of the gas or 
stellar emission. We obtained i by matching the distribution 
of [O III] and Ha emission with an ellipse. From this ellipse 
we also got the PA of the emission. We corrected the veloc- 
ity maps for the systemic velocity, using the values obtained 
in the central spectrum, and finally we could run ROTCUR. 
We fixed the following input parameters: inclination (i = 
40°), which was the same for the [O ill] and the Ha maps, 
the systemic velocity (0.0 km s~^), the expansion velocity 
(0.0 km s~^) and the coordinates (x, y) of the center of the 
galaxy. The other input parameters were the radii of the 
concentric rings in arcsec, the first guess for the rotational 
velocity of each ring and the position angle, and we let them 
free. To begin, we applied this task to the whole 2D velocity 
map, fitting both the approaching and the receding side at 
the same time. This output was used as input for the task 
VELFI which reconstructs the model of the observed veloc- 
ity field, then we checked if the model could well reproduce 
the velocity map. On the other hand, to get more accurate 
values of the deprojected velocity in order to build the ve- 
locity curve, the velocity maps were fitted in two steps: first 
the receding and then the approaching side. We compared 
the gas rotation curves with the stellar rotation curve, ob- 
tained applying ROTCUR to the stellar velocity map (i = 54° , 
PA^45°). We found that the kinematical behaviour and the 
velocity values are similar for stars and gas. This implies 
that the kinematics of the ENLR gas and of the stars are 
dominated by the same gravitational potential (see Fig. [T5|) . 
However, looking at the 2D maps we have already seen that 
the kinematical axis of stars is inclined of about 30° with re- 
spect to that of the gas fFig. [T3|) . This could be a projection 
effect: gas and stars could be distributed in a different way 
and in a different plane. Stars are in agreement with the 
continuum emission, while gas velocity maps are oriented 
according to the gas emission. We estimated the errors in 
velocity by measuring the position of some nightsky lines 
and assuming that the main contribution to the errors is 



200 




-200 - 



Figure 15. Rotation curves for gas and stars. The solid line is 
the best-fitting of the gaseous component. The median errors are 
12 km s"-*^ for [O ill] and 19 km s"-*^ for Ho;. The spatial resolution 
is 1 arcsec = 516 pc. 



due to wavelength calibration errors. We found median val- 
ues of 12 km s~^ for [O III] and 19 km s~^ for Ha. No error 
is reported in the output of starlight for stellar kinemat- 
ics. In principle, we s hould take into acco unt the S/N ratio 
of the emission lines (|Corsini et al.l ll999'): when the S/N is 
high, the error is due only to errors in calibration, instead 
when the S/N is low, the error is due to the badly fitted 
line profile. In our case, the S/N of these emission lines is 
high in the whole field. The errors may be larger in presence 
of asymmetries. For the high resolution data we also have 
the error due to the fitting of the two components. By as- 
suming that gas and stars are on circular orbits in a plane, 
we c an fit the rotation curves with the following formula 
(Bert ola et"al| ll991): 



Vc{r) 



(1) 



(r2 + cg)p/2 

where A, co, and p are parameters and r is the radius in 
arcsec. The p parameter is between 1 and 3/2. We fit the gas 
rotation curve, using both [O III] and Ha (see Fig. [15} with 
A = 300 km s~^, co = 2.5 arcsec and p = 1.1. The fitting of 
the stellar and gas rotation curve show a good agreement. 
By assuming a spherical potential, we can infer the mass 
distribution M(r) using the Virial Theorem. Considering as 
maximum extension r — 7 arcsec (corresponding to 3.6 kpc). 
Within this radius, we found M = 4.45 x 10^° M©. 

In order to perform a more accurate analysis of the kine- 
matical behaviour of the ionized gas, we focused on the [O ill] 
A5007 line of the higher resolution MPFS spectra. By using 
a multiple Gaussian fitting, we identified two kinematical 
components, a narrow and a broad one, which suggest a 
complicated kinematic structure of the NGC 7212 ENLR. 
We decided to apply the mid AS package XAlice, which is 
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Figure 16. Results of the multi-Gaussian fitting. Top: [O ill] ve- 
locity map obtained for the the main component (left) and the 
broader one (right). Bottom: FWHM maps for the main compo- 
nent (left) and the broader one (right). 



more effective than PAN in this case, and we obtained the 
flux, velocity and FWHM maps for both narrow and broad 
components (see Fig. [16]). The two components have dif- 
ferent velocities and different inclinations of the minor axis. 
For the narrow component we found values ranging between 
-120 km s"^ and +170 km s"^ and FWHM values up to 
350-400 km s~^. The broader component shows lower veloc- 
ity values, from -80 km s"^ to +120 km s"^ and FWHM 
up to 700-750 km s"\ We estimated PA = 130° for the 
stellar map, PA = 97° for the narrow [O ill] component and 
PA = 140° for the broad one. However, they are not aligned 
with the stellar velocity, as well. The high turbulence re- 
gion observed in the single-fitting [O ill] map and oriented 
as the minor axis of the gas emission is still present in both 
components. 

Having found hints of multiple components in [O ill] 
profiles, we decided to explore the gas kinematics, by 
analysing a high resolution echelle spectrum (R^SOOO, in- 
strumental FWHM'^0.85 A) obtained with the slit oriented 
along the ionization cones. By simply looking at the spec- 
trum it is clear that the emission is due to various sub- 
components at different velocities. We can see from Fig. [T71 
that the emission line profile of [O ill] is strongly variable 
and depending on the spatial position. Notwithstanding high 
spectral resolution, the components of the emission line pro- 
files of the high ionization gas are still broad, therefore we 
decided to analyse the lower ionization gas, using the [Sii] 
doublet, in which the emission lines are narrower and it is 
easier to identify and distinguish the different components. 
We found at least four components at different velocities (see 

Fig. [HI): vii 300 km s"\ vi2 - -150 km s"\ vis = 

km s"^ and vi4 ^ 170 km s"\ for the [Sii] A 6716. These 
components correspond to the ^721, V22, V23 and ^724 of the 
[S II] A6731. The FWHM of the V13 component of [S 11] A6716 
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Figure 17. The [Olll] emission line profiles at different distances 
(in arcsec) from the nucleus, extracted from the echelle spectrum. 

is about 200 km s~^. Maybe the lower velocity components 
are blended in the broad [O ill] emission line profiles and are 
not detectable in lower resolution data. We found an higher 
velocity component for the [S 11] A6731 at V25 450 km s~^, 
which is not visible in the [S 11] A6716 profile, because it is 
weak and it falls inside the [S 11] A6731 profile. It is worth 
noting that the ratio between the emission lines components 
changes with the distance from the nucleus, suggesting dif- 
ferent physical conditions of the ionized gas. In fact, in the 
HST image of NGC 7212 it is clear that the ENLR is made 
of filaments, likely regions at different velocities, which are 
probably characterized by different physical parameters. 



5 IONIZATION CONE 

From the ionization map we have detected an elongated and 
extended highly ionized structure with high values (10-14) 
of [Oiii]/H/3, up to 4 kpc from the nucleus, pointing out 
for the first time the presence of extended ionization cones 
in NGC 7212 (see Fig. I19|) . In our map, this structure is 
clearly visible both North- West and South - East of the nu- 
cleus, while in iFalcke, Wilson, Simpsonl (|l998l l only the 
south-eastern side is detected, probably because of the dust 
lane well visible in our B — R image. In fact, they calcu- 
lated the [Qui] /Ha ratio withou t any reddenin g corre ction. 
The same side was detected by ISchmitt et aD (|2003l ) with 
[O III] HST images. The lower limit of their measured flux 
is ^ 6.5 X 10~^^ erg s~^ cm~^ arcsec"^, while our de- 
tection limit is about 2 x 10~^^ erg s~^ cm~^ arcsec"^. 
burret k WarinI (|l990l ) measured [O iii]4959+5007/H/3 -19 
in the nucleus of NGC 7212. In their spectra this ratio 
ranges from 5 to 28 . The reddening corrected value from 
iBennert et al.l (|2006l l is — 16 in the nucleus and varies be- 
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Figure 19. The ionization map with the cone contours overlaid (left), the WFPC2 image of NGC 7212 (middle) and the same with the 
cone contours overlaid (right). 




Figure 18. [Sll] emission line profiles at different distances from 
the nucleus. The peaks corresponding to the various components 
are consistent. 



Figure 20. The [Olll]/H/3 ratio as a function of the distance 
from the nucleus. The squares (triangles) are the spaxels inside 
(outside) the cone. 



tween 6 and 17 in the central regions. Taking into account 
the [Oiii] A4959 contribution, we found [Oiii]/H/3 = 16 in 
the nucleus, values ranging between 12 and 19 at about 
7 arcsec North- West of the nucleus and between 12 and 
16 at about 5.5 arcsec South-East of the nucleus. We se- 
lected the spaxels of the cone considering only those show- 
ing [Oiii]/H/3 > 8. The aperture angle of the cone is ^70°. 
As we have already mentioned, this angle could be larger 
to account for the regions located far from the nucleus in 
a direction orthogonal to the cones, but showing AGN-like 
ionization. The cone has a total size of 12 arcsec correspond- 
ing to ^-^6 kpc and it is oriented at PA = 150°. 

The FWHM of the ionized gas inside the cone is around 
200-300 km s"^ The density is < 10^ cm~^, without a ra - 
dial trend, contrary to what found bv IPadina et al ] (|2010l ) 
in NGC 5252. These authors derived a relation from 
the constant [Oiii]/(0.5-2 keV) flux ratio, which suggests 
a radial-independence of the ionization parameter. We ob- 



served a gradient of temperature, with higher values in the 
northern region. We plotted the [Oiii]/H/3 ratio versus the 
distance from the nucleus, dividing the spaxels inside and 
outside the cone (Fig. [20)). As expected, we found that out- 
side the cone the ratio decreases with the distance from the 
nucleus, while inside the cone the ratio is high even far from 
the nucleus (up to 4 kpc). In order to understand if the AGN 
can account for such a high level of ionization, we made 
some energy balance considerations. Firstly, we calculated 
the L(Ha) from the reddening corrected Ha flux. Then the 
number of ionizing photons needed to pr oduce such lum i- 
nosities, Npu = 7.3 x lO^^L(Ha) ph s"^ (lKennicuttlll998l ). 
was estimated and compared with the number of photons 
emitted by the active nucleus, diluted by the covering fac- 
tor Q/47r, where Q is the solid angle of a region as seen 
from the nucleus. We did not make any hypothesis about 
the spatial distribution of the gas within each region. By 
assuming that the nuclear source is an isotropic emitter of 
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radiation absorbed along the line of sight, we estimated the 
nuclear ionizing photons by averaging the calculated number 
of photons ionizing the regions surrounding the nucleus and 
located within the cones, after having removed the effect of 
the covering factor. We found N^^j^"" = 9.4 x 10^^ s-^ The 
median value of the ratio between the observed value and 
the diluted one is 0.84, indicating that the active nucleus is 
the dominant ionizing source of the regions within the cone. 

Estimating the mass of gas can give an idea of the origin 
of the gas, in fact a large amount of gas can be explained 
with the acquisition of external material resulting from a 
merger or an interaction event. The mass of the gas in the 
ionization cone was estimated following two methods. The 
firs t one makes use of the Ha luminosity, as explained by 
[h3 (|2009l ): 

the second one assumes a Galactic dust-t o-gas ratio and 
makes use of the interstellar extinction IFu &; StocktonI 
(|2007l V We used a modified version of their formula: 

M{Av) = 1.5 X 10^kpc"Vag"^Ay(s x O^Mq (3) 

where s is the scale in kpc arcsec"^ and is the size of a 
region in arcsec. The results of the two methods are dif- 
ferent: in particular, the second method gives larger val- 
ues. The total values are M{Lh(x) — 5.16 x 10^ M© and 
M{Av) = 3.16 X 10^ M©. The reasons of this discrepancy 
could be a combination of two effects. First, the A^e value 
used in Eq. [2]is generally based on the [S ii]6716/[S ii]6731 
ratio, which fails in case of low electron density. Second, 
the flux of recombination lines depends on electron density 
squared. Therefore, the first method can cause an underes- 
timate of the ENLR mass, which is likely constituted by a 
large fraction of low density gas. 

We measured the asymmetry of t he [O III] e missi on lines 
by applying the method published bv lWhittld (|l985l ) to the 
higher resolution spectra. We calculated the area under the 
emission line and the wavelength corresponding to the 10 
per cent (Aio), 50 per cent (A50) and 90 per cent (Ago) of 
the emission line area. Then we obtained the asymmetry 
asym= where a — A50 — Aio and b — Ago — A50. In case 
of a blue wing in the line profile, the asymmetry is positive, 
while when a red wing is present, the asymmetry is negative. 
We reconstructed the asymmetry map (Fig. [21)) which shows 
two distinct regions, one with positive (C region, in Fig. 
[2T]) . and the other with negative asymmetry (A region, blue 
colours in Fig. [2T)) . These two regions match the ionization 
cone contours. We found [O ill] emission line profiles with 
blue wings in the northern part of the ionization cone, while 
the southern part of the cone is characterized by [O ill] with 
red wings. In the nucleus we found a symmetric profile. This 
suggests the presence of gas in radial motions, outflow or 
inflow, inside the cone. The asymmetry is not correspondent 
to the higher FWHM regions, that is orthogonal to the cone, 
therefore these high values of FWHM are not an effect of 
multiple kinematical components. 
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Figure 21. Top: the asymmetry map with the cone contour over- 
laid. Bottom: three examples of [O ill] profiles in different regions 
of the FoV. Panel A shows a spectrum with asymmetry equals 
to —0.14, panel B shows a profile without asymmetry, panel C 
shows a spectrum with asymmetry=0.21. 



6 CONCLUSIONS 

We studied the physical and kinematical properties of the 
circumnuclear gas in the nearby Seyfert 2 galaxy NGC 7212. 
We analysed this object by means of integral field and echelle 
spectra and broad-band images. We pointed out for the first 
time the presence of an extended ionization cone in NGC 
7212, with high values of [Oiii]/H/3 (up to 12), at a large 
distance from the nucleus (up to 3.6 kpc). The cone is ori- 
ented NW-SE, at PA = 150° close to the photometric minor 
axis of the galaxy, with an opening angle of about 70° . The 
cone is more extended in the North (7 arcsec = 3.6 kpc) than 
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in the South (5.5 arcsec = 2.8 kpc) direction, while in HST 
[O III] pubhshed images only the southern high ionization 
emission was detected, showing a structure made of clouds 
or filaments. The evidence of dust located in the ENLR also 
supports this idea that the gas in the cone has a filamen- 
tary structure, and suggests that a potential source of error 
raises in the dereddening, because how the gas and dust are 
mixed in and between filaments is in fact unknown. The 
mass of the ENLR was calculated by means of two different 
methods based on the Ha luminosity and on the interstellar 
extinction, and it is likely between 5x10^ and 3x10^ M©. 
NGC 7212 is in an interacting triplet, with two galaxies in 
a clear on-going merger. We found high values of [Nil] /Ha 
and [Sii]/Ha towards the interaction region: suggesting a 
possible combination of ionization by the active nucleus and 
by shocks. From the 2D velocity maps, we found kinematical 
decoupling between the stars and gas. The velocity fields are 
misaligned, with the stellar kinematical minor axis aligned 
with the photometric minor axis, while the gas kinematical 
axis is tilted 30° with respect to the stellar one. We stud- 
ied the asymmetry of the emission line profiles inside the 
ionization cone and we found [O III] with blue wings in the 
northern side of the cone, while the southern side of the cone 
is characterized by [O ill] with red wings. In the nucleus the 
profiles are symmetric. This suggests the presence of gas in 
radial motions, which is confirmed also by the analysis of 
high resolution spectra. In fact, studying the echelle data, 
we found that the ionized gas is characterized by multiple 
kinematical components at different velocities. We found at 
least four components: at —300 km s~^, —150 km s~^, 170 
km s~^ and 450 km s~^, with respect to the recession veloc- 
ity. The ENLR gas metallicity was estimated by measuring 
the observed emission line ratios and comparing them with 
CLOUDY models, obtaining indications of sub-solar values. 
All these properties support the idea of an external origin 
of the ENLR gas in NGC 7212, likely due to gravitational 
interaction effects in act in this triple system. 
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